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Abstract. A promising practical application of entanglement is metrology, where quantum states 
can be used to make measurements beyond the shot noise hmit. Here we consider how metrology 
schemes could be realised using atomic Bose-Einstein condensates (BECs) trapped in different 
potentials. In particular, we show that if a trapped EEC is rotated at just the right frequency, it can 
undergo a quantum phase transition characterised by large-scale entanglement spreading across the 
system. This simple process of stirring can generate interesting quantum states such as macroscopic 
superpositions of all the atoms flowing in opposite directions aroimd a ring-shaped potential. We 
consider different trapping potentials and show how this leads to different entangled states. In 
particular, we find that by reducing the dimensionality of the system to one or two dimensions, 
it is possible to generate entangled states that are remarkably robust to the loss of atoms and so are 
ideally suited to precision measurement schemes. 
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1. INTRODUCTION 

One of the most exciting recent applications of quantum information theory has been to 
the field of metrology. A number of theoretical proposals [1, 2, 3, 4] have demonstrated 
that it is possible to use entanglement to surpass the so-called shot-noise limit. This limit 
is due to unavoidable statistical fluctuations that arise when discrete independent events 
are used to make a measurement. It is characterized by a precision that scales as 1 /\/N 
where N is the number of discrete events recorded, such as tosses of a coin or particles 
detected. 

If, instead, we remove the constraint that the subsystems are independent of one 
another and allow them to be entangled, then it is possible to reach the Heisenberg 
limit whereby the measurement precision scales as 1/A'^. There have even been some 
proposals for schemes that exploit nonlinearities to go beyond the Heisenberg limit 
[5,6]. However, it is now thought that these schemes misinterpret the resources required 
and so the Heisenberg Umit is indeed fundamental [7]. 

Metrology schemes are often based on interferometers because they offer unrivalled 
precision in the measurement of phase shifts or path length differences. As well as their 
familiar linear forms, interferometers can also be successfully implemented in ring ge- 
ometries to make accurate measurements of angular momentum. These interferometric 
gyroscopes form a key component of many modem navigation and stabilisation systems. 
They work by exploiting the different path lengths experienced by light as it propagates 
in opposite directions around a rotating ring. For instance, in the Sagnac geometry [8] 



photons are put into a superposition of travelling in opposite directions around a ring 
and, when rotated, the two directions acquire different phases due to their different path 
lengths. This phase difference is directly related to the rate of rotation and can be mea- 
sured by recombining the two components at a beam splitter. 

The precision of gyroscopic devices can be improved further still with the use of 
entangled atoms rather than photons. The precision with which a Sagnac interferometer 
can measure a phase difference is [9] 

4n 

A0 = —CIA, (1) 
Av 

where A is the particle wavelength, v is the particle speed, Q. is the angular velocity of the 
ring, and A is the enclosed area. For equal values of D. and A, the ratio of the precision 
for atoms with mass M to photons with frequency CO is Mc^/h(D, which for reasonable 
parameters is about 10^ ^ So, even though the particle fluxes are likely to be much lower 
for atoms and the enclosed areas smaller, the very large enhancement factor for atoms 
means that they are likely to offer considerable advantages over photons. 

The use of entangled atomic states to make precision measurements of rotations has 
received a lot of recent attention. It is possible to achieve sub-shot-noise measurements 
by using number- squeezed atomic states as the input to an interferometer and so there 
has been a lot of research into the generation and uses of these states [10, 11, 12, 13]. 
Several proposals have already been made to use different non-classical atomic states 
to make general phase measurements with sub-shot noise sensitivities [14, 15, 16, 
17, 18, 19, 20]. It has also been shown that uncorrelated atoms can also achieve sub 
shot-noise sensitivities of rotational phase shifts [21] by using a chain of matter wave 
interferometers, or a chain of gyroscopes. 

It is clear that there is a link between entanglement or squeezing and the measurement 
precision that can be achieved, however the precise relationship is not obvious. While it 
can be shown that a NOON state of the form, 

\llf) = ^ma\0)b+\0)a\N)b), (2) 

where the two kets represent Fock states of two states a and b respectively, gives the 
best possible precision (i.e. it saturates the Heisenberg limit) in the lossless case, the 
same is not true when we account for the loss that will inevitably be present in any 
realistic system. The problem is that NOON states are known to be extremely fragile. 
The loss of a single particle destroys the superposition. So while, they have excellent 
short-term sensitivity, their long-term stability is very poor. This means that they can 
only be used to make very short timescale comparisons, which degrades their overall 
precision. Think of frequency standards: the longer we can compare our system with the 
standard, the more accurately we can detect any difference in their frequencies. So there 
is an important trade-off between the entanglement in a system and the robustness of the 
state to decoherence. 

We can quantify how accurately a state can measure an unknown parameter with 
the quantum Fisher information and its associated Cramer-Rao bound. Quantum Fisher 



information is defined as, 

Fg = Tr[p(0)A2], (3) 

where p(0) is the density matrix of the system, ^ is the parameter to be measured, and 
A is the symmetric logarithmic derivative defined by, 

^ = ^[AP(0)+P(0)A]. (4) 

In the eigenbasis of p(0) this is, 

where Xij are the eigenvalues of p(0) and p'{^) — dp{^)/d^. If Xi + Xj — then 
{A)ij — 0. The Fisher information is related to the uncertainty of the measurement via 
the Cramer-Rao bound, A0 > 1/y^^. For a pure state, | V^(0)), we have, 

fq = 4 [(v/(0)|i/(0)) -\{xi/mwm?] , (6) 

where = d\\i/{(j))) /d^. Fisher information will serve as our yardstick for com- 

paring the performances of different entangled states in metrology. 

In this paper we will discuss how a variety of entangled states that are useful for 
metrology can be created simply by stirring Bose-Einstein condensates (BECs) trapped 
in potentials with different geometries. We start with the case of NOON states and then 
show how, by modifying the potential and dimensionality of the system, it is possible 
to create different entangled states that perform better in the presence of decoherence. 
The general study of the relationship between entanglement, phase transitions, and 
the geometry of the trapping potential is an exciting new area of research that is just 
starting to open up. This development has been motivated, in part, by recent outstanding 
experimental advances in the ability to trap atoms in a wide variety of time-varying 
potentials. 



2. LATTICE AND CONTINUOUS RING POTENTIALS 

The first trapping geometry that we will consider is a lattice of potential nodes in a 
ring configuration. For simplicity, we will consider a lattice consisting of just three sites 
as shown in Figure 1(a) since this is the minimum number required for a ring. The 
different sites are coupled to one another via quantum mechanical tunnelling through 
the barriers separating them. The atoms on any given site also experience a nonlinearity 
due to collisions with other atoms on that site. This trap geometry has already been 
experimentally demonstrated [22, 23]. 

Assuming that the system is sufficiently cold that we can consider only one trap 
level at each site, the system can be described by the Bose-Hubbard Hamiltonian. The 
rotation of the sites can be incorporated into the Bose-Hubbard Hamiltonian [24] simply 




FIGURE 1. (a) The three site ring lattice. Atoms are trapped at sites a, b, and c and can tunnel between 
adjoining sites, (b) Energy levels (in units of J) plotted as a function of the rotation, (j),foi J = U =1 and 
= 3. A clear anti-crossing between the two lowest energy levels is seen for (j) = n. 



by including appropriate phase factors in the coupling terms. This gives the 'twisted' 
Hamiltonian, 



H 



-J 



a^b + b^c + c^a ) -\-h.c. 



+ U[a'^a^ + b^^b^ + c^^c^], (7) 



where a, b and c are the annihilation operators of atoms in the three sites, U is the on-site 
interatomic interaction strength and J is the coupling between adjacent lattice sites due 
to tunnelling through the barrier. The phase factors in the coupling terms are known as 
Peierls phase factors. 

It is convenient to rewrite the Hamiltonian in the quasi- momentum basis {a, j8 , 7} 



a 



1 

7! 



{a^b-\-c), 



/3 = ^(a + Z7e'2^/3+ce'4'^/3), 
V3 



(8) 



These new basis states respectively correspond to zero flow, one quantum of clockwise 
flow, and one quantum of anticlockwise flow, where throughout this paper we will use 
the convention that a positive phase variation corresponds to clockwise flow. 
In this basis, the twisted Hamiltonian has the form 

H = -7{(2a^a-/3'^/3-7V)cos(<|)/3) + V3(/3^/3-7V)sin(0/3)} 
U 



(9) 



The full energy spectrum of the system can be found by diagonalising this Hamiltonian. 
However, we can gain a lot of insight about the behaviour of the system simply from 
inspection. When there are no interactions {U = 0) the second and third line of the 
Hamiltonian vanish and when there is no rotation ((^ = 0) we see that the ground state is 
the non-rotating state, a, and the other two states are degenerate, as we would expect. If 
we were to increase the rotation rate to = 2tz, the ground state would be all the atoms 
in /3, i.e. all rotating with one quantum of angular momentum. In between these values, 
i.e. when ^ = K, a and /3 become degenerate ground states. However, if we include 
interactions, we see from (9) that the different modes are now coupled. This serves to 
lift the degeneracy and causes an anticrossing at ^ = :;r. This can be seen in Figure 1(b), 
where the energy levels are plotted as a function of for J = U = I and three atoms. 

To a very good approximation, the ground state at the critical rotation rate = n) 
is a symmetric superposition of all the atoms in a and all in /3, i.e. a macroscopic 
superposition or NOON state, 

\¥) = ^{\N)a\0)p + \0)a\N)p). (10) 

This is an intriguing result because it means that, simply by stirring a Bose condensate 
at just the right rate, it is possible to generate a highly entangled state. 

Strictly speaking, when the barrier heights are equal, we only get a NOON state by 
this method when we have commensurate filling, i.e. the total number of atoms is an 
integer multiple of the number of lattice site. However, it has been shown that we can 
get around this problem simply by ensuring that the barrier heights are slightly different 
[25]. 

Very similar results have also been demonstrated for atoms trapped in a continuous 
loop with a small barrier at one point [26] as shown in Figure 2(a). The atoms can tunnel 




FIGURE 2. (a) The continuous loop potential with a small barrier that the atoms can tunnel through, 
(b) Energy levels (in units of Ti^/ (2Mr^), where M is the mass of an atom and r is the radius of the loop) 
plotted as a function of the rotation, (j). The dashed lines are for the case without the barrier and the solid 
lines are the energy levels when the barrier is included. An anti-crossing between the two lowest energy 
levels is seen for (j) — 7Z, i.e. when the loop is rotated at a rate corresponding to half a quantum of angular 
momentum per particle. 



through the barrier and the loop is rotated. The energy levels for this system with and 
without the barrier present are shown in Figure 2(b). At a rotation rate corresponding 
to a phase ^ = 7t around the loop (i.e. when the loop is rotated at a rate corresponding 
to half a quantum of angular momentum per atom) there is an anti-crossing for the case 
that the barrier is present. As in the case of the lattice, the ground state at this critical 
frequency is a macroscopic superposition of all the atoms not rotating and all rotating 
with one unit of angular momentum. 

Both the lattice and continuous loop geometries discussed above could be used as 
entanglement enhanced gyroscopes. The general idea is as follows [27, 28]. We begin 
with a BEC trapped in a non-rotating version of one of these potentials. We then 
abiabatically increase the rate of rotation until we reach the critical frequency. At this 
point, the system should be in the ground state, which is a NOON state of rotating and 
non-rotating components. Next we make a sudden (non-adiabatic) change to the rotation 
frequency and allow the system to evolve for a while. Finally we adiabatically reduce the 
rotation rate well below the critical frequency and detect the atoms in the rotating and 
non-rotating basis. In the ideal case, all atoms should be found in one or the other. The 
probability of these two outcomes depends on the evolution after the non-adiabatic shift 
and should provide a measure of the angular momentum of the system. Importantly, the 
rate at which these probabilities change with the angular momentum should be A'^ times 
faster than if we used the usual Sagnac setup where the atoms are not entangled, but each 
is put into a superposition of going each way round the loop. This enhancement of the 
frequency of the interference fringes due to entanglement has already been observed in 
optical systems [29, 30]. It offers the possibility of making measurements with enhanced 
sensitivity right up to the Heisenberg Umit. 



3. DECOHERENCE 

So far we have considered only NOON states. The Fisher information for these states 
given by (6) is Fq — N^, which means that they can reach the Heisenberg limit, A<^ > 
1 /N when there is no decoherence. However they are known to be very fragile and so 
quickly lose their advantage when any loss is present. For more than a small amount of 
loss, the case of unentangled atoms each put into a superposition of the modes a and j8, 
outperforms the NOON state. 

However, it is possible to find states that embody the best aspects of each of these 
cases, i.e. states that outperform the precision of unentangled particles and yet don't 
suffer from the same fragility as NOON states. One such example, known as a bat state, 
is formed by passing a dual Fock state, | — \N) \N) through a 50:50 beam splitter 
[2, 31]. It has the form | i/a) = Zm=o^m\'^m)\2N -2m) where, 

_ V(2m)!(2A^-2m)! 

2^m\{N-m)\ ' ^ ^ 



These coefficients are plotted in Figure 3(a). We see that the state is quite like a NOON 
state in that it has a large number variance. This means that it has a large Fisher 
information and hence can measure phase shifts with high precision. However, unlike a 




FIGURE 3. (a) Plot of the coefficients of a bat state with 10 atoms, (b) Comparison of the quantum 
Fisher information for three different states (each initially with 10 atoms) as a function of the fraction of 
particles lost: (i) NOON state, (ii) the bat state shown in (a), and (iii) the case where the atoms are not 
entangled with one another and each is in an equal superposition of modes a and j3 . 



NOON State, if we detect which mode one of the atoms is in, this does not betray which 
mode all the remaining atoms are in. In other words, the bat state should be much more 
robust to the loss of atoms than the NOON state. 

In Figure 3(b) we have plotted the Fisher information for these three different states 
initially containing 10 atoms as a function of the fraction of particles lost. For simplicity, 
we have assumed that there is equal loss for both of the modes. We see that in the lossless 
case, the NOON state has the highest Fisher information, as expected. However, for 
moderate loss, the bat state outperforms both the NOON and unentangled states. These 
differences become even more dramatic as the number of particles is increased. 

Since loss is inevitable in any realistic system, it would seem that it would be very 
useful if we could create bat states with rotating matter waves. A scheme to do this 
was recently proposed in [28]. However due the highly squeezed state required to create 
the bat state and the unwanted effects of particle interactions in the scheme, this matter 
wave gyroscope is limited to a small number of particles thereby limiting its practical 
precision capabilities. So the question is, can we find practical ways for creating other 
entangled states that are well-suited to metrology in the presence of loss? This is still 
very much a new field of enquiry, however two intriguing examples show how this is not 
only possible, but that we can do even better than bat states. This is achieved by reducing 
the trapping potential down to one or two dimensions. In the following two sections we 
review these two cases. 

4. TWO-DIMENSIONAL PANCAKE TRAP 

Dagnino et al. recently demonstrated a very interesting model for creating entangled 
atoms in a rotating potential [32]. They considered a two-dimensional 'pancake' shaped 
potential. One of the dimensions can be removed in practice by taking the trapping 
frequency in (say) the z-direction to be much greater than trapping frequencies in the x 



and y directions, i.e. (Oy ^ (Ox,y. If ho)^ is sufficiently large compared with the interaction 
energy, the dynamics along z are frozen and the gas is effectively two-dimensional at 
sufficiently low temperatures. 

The Hamiltonian for this system is. 



where the sums are over the particles, M is the mass of an atom, €l is the rate of rotation 
of the potential, L^i is the z-component of the angular momentum of particle z, and (Oxy 
is the radial frequency of the trap in the xy plane. The interactions between atoms are 
parameterized by the unitless quantity, g — a^%%MoiJh, where a is the 3D scattering 
length. The asymmetry is given by Yi = 2AMo)xy{x^ ~y^)^ where the coefficient A 
measures the strength of the asymmetry. 

When A ^ 1 and Q. ^ (Oxy, the single particle energy levels are grouped into Landau 
levels separated by energy + (Oxy). This means that if the interaction energy is small 
compared with this spacing, to a good approximation, the system is restricted to the 
lowest Landau level [33]. A basis for the single particle states of the lowest Landau level 
for A = is. 



¥m=J^(^y (x + /3;^e-^--(-^+^^)/2^ (13) 
y rniTt \ h J 

where m > is an integer. Using this basis, we can rewrite the Hamiltonian (12) in 
second quantised form. 



Pico i4. / 

H = hcOxyN + h{(Oxy-i^)L-\ ^ ^ ^ ^/(m + 1 ) (m + 2)aj„am+2 + \/m{m-2)al^a„ 



h(Oxyg Y Y 5mi+m2,ni+W2 {m.i+m2)\ 

A-n- 1^ 1^ 



t t 



4n «^^,„f^,Vmi!m2!ni!n2! 2'«i+'«2 

where aj is the annihilation operator for a particle in the state Xj/j given by (13), and 
N = Y^m^tflm and L = Y^m^^^rflm are respectively the number operator and the total 
z-component angular momentum operator. 

The spectrum of energy levels for different parameters can be found by directly 
diagonalising (14). In the case that the potential is symmetric (i.e. A = 0), the two lowest 
energy levels become degenerate at some critical rate of rotation, Q.c=\ — gN / If 
we include some small asymmetry, this degeneracy is lifted and there is an anti-crossing 
at the critical frequency. This should be compared with the two other geometries we 
have discussed: the lattice loop and the continuous loop. In the former, the degeneracy 
was lifted by including interactions between the atoms and in the latter it was lifted by 
including a small barrier into the loop. 
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FIGURE 4. The ground state coefficients of the asymmetric 2D potential at the critical rotation rate, 
where the state has the form of Eq. (15). The parameters chosen for this plot ais N = 12, g = 0.5 and 
A ==0.03. 

At the critical point, the ground state of the system (for even N) can be written 
approximately as 

N/2 

\xir)=Y^Cnr\2m)\N-2m), (15) 

111=0 

where the first ket represents an even parity state that is approximately a coherent 
superposition of the states i/Zq and 1//2, while the second ket is an odd parity state well- 
approximated by yAj [32]. The coefficients of this ground state at the critical rotation rate 
are shown in Figure 4 for the parameters = 12, ^ = 0.5 and A = 0.03. We see that the 
state is qualitatively similar to the bat state shown in Figure 3(a). In particular, it has a 
large number variance and so has a large Fisher information meaning that it should be 
able to detect phase shifts very accurately. Also, the fact that the detection of a particle in 
one of the two particular modes that make up the state does not destroy the superposition 
means that it should be relatively robust to loss. 

Overall, this seems to be a relatively straightforward scheme for generating an en- 
tangled state that is well-suited to metrology in the presence of loss. Developing a full 
scheme for how this system could be exploited as a quantum-limited gyroscope as well 
as analysing the behaviour of the state in the presence of decoherence are interesting 
directions for future work. 

5. ONE-DIMENSIONAL LOOP AND THE TONKS-GIRARDEAU 

LIMIT 

Another promising example occurs in a one dimension ring trap [34, 35]. Suppose we 
have a ring potential with a moving barrier, similar to that shown in Figure 2(a), but 
where the potentials are so tight in the radial direction that we are left with only the 
angular, 6, dimension. If the interactions between the atoms are strong, we enter the 
so-called Tonks-Girardeau (TG) regime where the atoms cannot pass one another. 



The Hamiltonian that describes this system is 



N 
i=l 



[ .d ay ^ 



(16) 



where the sum is over all atoms of mass M. The circumference of the ring is L, 
which gives x — 6L/2k as the atom's position on the circumference of the loop and g 
is the effective one-dimensional interaction strength between the atoms. The smallest 
nonzero kinetic energy of a single particle Eq = In^h^ / (ML^) provides a natural unit of 
energy for this system. The barrier is assumed to be narrow, so it can be described by 
a 5 -function and has strength b. The barrier stirs the atoms in the ring with tangential 
velocity v = Ml/ (ML) along the circumference of the ring and the Hamiltonian (16) is 
formulated in the co-rotating frame of reference. 

It is convenient to write the Hamiltonian in the second quantized form using angular 
momentum operators, 

'i'ix) = ^Y^e'^^'^/%, (17) 

where a| and create and destroy an atom with angular momentum kh, respectively. In 
this representation, the Hamiltonian is, 

H = Y^^oU-^^ «K + ^E4i^*2 + ^ E 4C^^i-A2+?' (18) 

k ^ ' ki,k2 ki,k2,q 

where the sums are over all quantized angular momentum modes, kh. Again the ground 
state and energy spectrum can be found by diagonalising this Hamiltonian, however, the 
Hilbert space is infinite so a truncated momentum basis has to be used. This means that 
the interaction strength needs to be rescaled as discussed in [34] . 

Much can be understood about the system by considering the terms in the Hamilto- 
nian. The kinetic energy term and the interaction term both conserve the total angular 
momentum and it is only the barrier term that couples states of different total angular 
momentum. When the barrier is rotated at a rate D. = 7U,the energy of the and h angular 
momentum modes become degenerate, as described in section 2. In the non-interacting 
case this means states with to Nh total angular momentum become degenerate. When 
interactions are introduced the degeneracy is lifted except the and Nh states and at a 
certain interaction strength, and with a finite barrier height, the ground state becomes 
similar to a NOON state. Increase the interactions further and the NOON state is lost, 
however the ground state is still a superposition of to Nh total angular momentum, be- 
cause the interaction does not couple states of different total angular momentum. What 
we find is the states above this certain interaction still produce measurements with the 
same precision as the NOON state, however, their robustness is far better. 

In the strongly interacting TG regime the state of the system looks very different 
from a NOON state. The single particle momentum distribution for the and Nh total 
angular momentum states are spread out and overlap each other. This means that if an 




Fraction lost 

FIGURE 5. (i) The quantum Fisher information of a Tonks-Girardeau gas as a function of the fraction 
of atoms lost. This is compared with (ii) the NOON state and (iii) the unentangled state. In each case the 
state contains 5 atoms before loss. 



atom is detected it cannot be discerned which of the total angular momentum states it 
came from and so it does not destroy the superposition. The performance of the TG 
superposition in the presence of loss is compared with the NOON and unentangled 
states in Figure 5 for A/^ = 5. The TG plot is for b/L — O-OOS^o (i-e. a small barrier) and 
strong interactions, g/L= 1085 / {2Tt)EQ, i.e. deep in the TG regime. The calculation was 
performed on a truncated basis with 18 angular momentum modes and the results do not 
change noticeably as the basis is increased. We see that the initial Fisher information is 
equal to that of the NOON state in the case without loss. When loss is considered, the 
reduction of the Fisher information is far slower and is always better than unentangled 
atoms. Another advantage of the system is the practical creation of the superposition. 
For the NOON state it was shown that even for the case of no decoherence it would be 
experimentally difficult to create for more than a few atoms due to the extreme precision 
required in the rotation rate [26] . The TG superposition is far more stable to experimental 
imperfections and energy level spacings between the ground and excited states are far 
larger allowing more rapid creation of the superposition. 

6. CONCLUSIONS 

The idea of creating different entangled states by rotating matter waves in trapping 
potentials with different geometries is an exciting new area of research. There are 
already some interesting results for how entangled states that are suited to quantum 
metrology schemes could be generated relatively easily. There is still a lot to be done. It 
will be particularly interesting to see how different symmetry properties in the potential 
manifest themselves in different features of the entangled states. This may lead to quite 
general quantum state engineering proposals for rotating matter waves. The models 



proposed here and their extensions may also prove to be valuable models for exploring 
the deeper connections between phase transitions, entanglement, and metrology. 
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